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Abstract. We have pursued two different methods to analyze the old stellar population near the Galactic plane, 
using data from the 2MASS survey. The first method is based on the isolation of the red clump giant population 
in the color-magnitude diagrams and the inversion of its star counts to obtain directly the density distribution 
along the line of sight. The second method fits the parameters of a disc model to the star counts in 820 regions. 
Results from both independent methods are consistent with each other. The qualitative conclusions are that the 
disc is well fitted by an exponential distribution in both the galactocentric distance and height. There is not 
an abrupt cut-off in the stellar disc (at least within R < 15 kpc). There is a strong flare (i.e. an increase of 
scale-height towards the outer Galaxy) which begins well inside the solar circle, and hence there is a decrease of 
the scale-height towards the inner Galaxy. Another notable feature is the existence of a warp in the old stellar 
population whose amplitude is coincident with the amplitude of the gas warp. 

It is shown for low latitude stars (mean height: \z\ ~ 300 pc) in the outer disc (galactocentric radius R > 6 kpc) 
that: the scale-height in the solar circle is /i z (i?Q) = 3.6 x W~ 2 Rq, the scale-length of the surface density is hn = 
0.42_Rq and the scale-length of the space density in the plane (i.e. including the effect of the flare) is H — 0.25-Rq. 

The variation of the scale-height due to the flare follows roughly a law h z (R) ~ hz(R@) exp ( [ 12 _ &R(kpc)] kpc ) 

(for R <~ 15 kpc; Rq = 7.9 kpc). The warp moves the mean position of the disc to a height z w — 1.2 x 
10~ 3 i?(kpc) 5 ' 25 sin(</> + 5°) pc (for R <~ 13 kpc; R Q = 7.9 kpc). 
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1. Introduction 

Star counts have been widely used for the study of 
Galactic structure (see Paul 1993), and are becoming in- 
creasingly important with the appearance of wide area 
surveys in the last decades (see for example Bahcall 1986, 
Majewski 1993, Reid 1993, Garzon et al. 1993, Price 1988). 
The use of progressively larger and more sensitive samples 
combined with detailed models of stellar galactic distribu- 
tion (for example Bahcall & Soneira, 1980, Wainscoat et 
al. 1992) has overcome most of the original uncertainties 
and difficulties in the analysis and interpretation of star 
count data. The advent of NIR detector arrays has permit- 
ted the detailed exploration of the stellar structure in hith- 
erto hidden areas of the Milky Way, such as the Galactic 
Plane and Bulge (Eaton et al. 1984, Habing 1988, Garzon 
et al. 1993, Hammersley et al. 1994, Ruphy et al. 1996, 
Lopez-Corredoira et al. 2000, Epchtein 1997, Skrutskie et 
al. 1997). This is because at these wavelengths interstellar 
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extinction is significantly less when compared to the vis- 
ible, while keeping the individual stellar contribution to 
the observed flux. 

There are, however, still controversial or totally un- 
known parameters in the description of the detailed stel- 
lar structure. Some of these are concerned with the radial 
and vertical distribution of the Galactic disc, and its spe- 
cific morphology. Several radial scale-lengths have been 
determined by different groups, even when using the same 
data sets. As a general trend, the shorter scale-lengths 
are found for the surveys using longer wavelengths (Kent 
et al. 1991), although it could equally be that different 
source type or regions of the Galaxy are being examined. 
There is far less debate on the vertical scale-height, with 
the canonical value for the old population in the Solar 
vicinity being about 300 pc. Still more work is needed, 
however, to differentiate between the distribution being 
an exponential or the sech 2 (van der Kruit 1988). Finally, 
features in the stellar disc like internal (Freudenreich 
1998, Lopez-Corredoira et al. 2001) or external (Habing 
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1988, Ruphy et al. 1996) cut-offs radii, flares (Kent et 
al. 1991, Alard 2000), warps (Djorgovski & Sosin 1989, 
Freudenreich 1998, Porcel & Battaner 1995, Alard 2000, 
Drimmel & Spergel 2001) and/or local corrugations of the 
galactic plane (Hammersley et al. 1995) are less well stud- 
ied, and are more controversial. 

In this paper we will address some of the above top- 
ics by analyzing the shape of the stellar Galactic disc in 
two ways. We will first extract a single stellar type by 
using the deep NIR color-magnitudes diagrams (CMD) 
from the 2MASS survey (Skrutskie et al. 1997). The anal- 
ysis will be based on the red clump giants, whose relatively 
high density and bright absolute magnitude give rise to an 
identifiable feature in the CMDs. Hence, the distance and 
extinction to each star can be directly determined only as- 
suming the absolute magnitude and color of the sources. 
As a second approach, we will use the star counts, also 
taken from 2MASS, and the results will then be com- 
pared with those obtained from the red clump method. 
Both methods will be used to examine the scale-length, 
scale-height, flare, outer disc cut-off and warp. 

The data for this work has been taken from the 2nd 
release of the 2MASS project (Skrutskie et al. 199 7, 



http: / /www.ipac.caltcch.cdu/2mass/releases/docs.htm| ). 
We have made use of the available data in that sec- 
ond release for near plane regions in the outer Galaxy 
(45° < I < 315°). 

2. Model of a non-warped Galactic disc 

In order to help to interpret the 2MASS star counts a 
simple model of the density distribution in the disc was 
built. 

The Galactic disc is assumed to be an exponential dis- 
tribution whose star density is: 



p(R, z) = p G exp ( - — ^ Q I < xp 



1*1 \ h z (Re) 



h z (R)J h z (R) 



(1) 



h s (R) = h z (R & )exp 



R — Rq 

h-RJlare 



(2) 



where p Q is the space density on the plane in the neigh- 
borhood of the Sun. 

Hence the density falls exponentially with the galac- 
tocentric radius R (scale- length hj{), and with the height 
z. We have preferred an exponential scale-height rather 
than a sech 2 , as Hammersley et al. (1999) shows that a 
sech 2 law did not give a good fit to the TMGS star counts 
near the plane. The effect of the flare is to increase the 
scale-height of the sources in the outer Galaxy, hence the 
scale-height (h z (R)) is dependent on R. We have assumed 
that the total number of stars at a specific galactocentric 
radius remains the same as it would have been without 
the flare and so all the flare is doing is to distribute the 
sources further from the plane. In the model this is in- 
cluded by having the scale-height increase exponentially 
with galactocentric radius and then have a normalizing 



factor ( h ^'^ j for the density. Whilst this is not the first 
model to include a flare, it is the first to use an exponential 
rather than a linear flare. 

Another way of expressing the density is: 



p(R, z) = p Q exp 



R — R(\ 



exp 



h z (R) 



where 
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(3) 



(4) 



hR. flare , 

i.e. H is the equivalent scale-length which is the combina- 
tion of the intrinsic scale-length of the disc with that of 
the flare. Written in this form the equivalent scale-length 
takes the place of the normalizing factor. The flare in- 
creases with R and so removes proportionally more stars 
from the plane with increasing distance. The net effect is 
that on the plane the star density falls far more rapidly 
with increasing galactocentric distance than the intrinsic 
scale-length predicts. Hence, hn is the scale-length of sur- 
face density, and H is the scale-length of the space density 
in the plane. 

We have assumed a distance from the Sun to the 
Galactic center to be Rq — 7.9 kpc (Lopez-Corredoira 
et al. 2000). 

Therefore, this model contains only three free parame- 
ters (ha, Ii z (Rq), ha fi are ) whose determination is carried 
out by minimizing the \ 2 °f the fit in the following way. 

For a given set of parameters, the space density on the 
plane in the neighborhood of the Sun can be determined 
using the formula 



E 



» a? 



(5) 



where di is the measured density in each region, cjj the 
poissonian errors in the counts and ij is the model predic- 
tion with pQ = 1. Hence, for each set of model parameters 
we calculate p@ through eqn. (^), and after that we cal- 
culate \ 2 for that particular fit. Once we have calculated 
all the x 2 f° r all se ^s in the parameter space, we select the 
one which gives a minimum x 2 ■ The errors are due to the 
errors in the counts (assumed to be poissonian) and errors 
coming from the calculation. 

The density equation can be further simplified on plane 
(b « 0), as long as \z\/h z {R) « 1 along the whole line of 
sight: 



p(R, z = 0) = p exp 



R — Rq 

H 



(6) 



If the external disc were truncated, then another pa- 
rameter that could be added is the radius of the cut-off 
-Rcut-off- It will be shown in the analysis presented here 
(§3.5), however, that there is no evidence for a cut-off to 



at least R — 15 kpc, and so this parameter will be unnec- 
essary in the fits. Beyond R « 15 kpc, there are very few 
stars detected and so very little can be said. 



M. Lopez-Corredoira, A. Cabrera-Lavers, F. Garzon, P. L. Hammersley: Disc/2MASS 



3 



We have ignored the contribution of the spiral arms 
since in the outer Galaxy their contribution to the star 
count is small and even then it will be at the brightest 
magnitudes being considered here, giving negligible con- 
tribution for ran < 14.0. 

Other possible additions to the standard disc model 
are large scale warps and the presence of the thick disc. 
The warp will be considered separately in §pl 

The thick disc (Gilmore & Reid 1983; Buser et al. 1998, 
1999; Ojha 2001) will not be studied here as a separate 
population since the analysis presented here is not sen- 
sitive to thick disc, as regions in general near the plane 
are used where its relative contribution is small. It will be 
shown in sections |3.5| , |] that the thick disc could consti- 
tute as much as a 10% of the total low latitude sources 
selected in this paper, so the derived scale-height will be 
somewhat larger than would be the case for the thin disc 
alone. Here we will use the term "old disc" as the aver- 
age structure traced by the old disc stars, including any 
thick disc component (mean height around 300 pc). The 
determined mean scale-height might be larger at higher 
latitudes if there is a thick disc. 



3. Using CMDs to obtain the parameters of the 
disc 

3.1. The selection of giants 

We have built (J - K, m K ) CMDs from the 2MASS data 
set for selected regions of the sky (Fig. |l|). We have taken 
the data 'as is', without further processing. The data have 
been grouped in regions cove ring between 2 and 5 squared 
degrees on the sky (see §3.4 for details). 

The position on a CMD of a star is determined from 
its absolute magnitude, intrinsic color, distance and ex- 
tinction. Any horizontal displacement is caused by extinc- 
tion alone, and hence if the intrinsic color of the source 
is known, the total extinction can be directly determined 
from the measured color. The vertical motion, however, is 
caused by both distance (doubling the distance makes the 
source 1.5 mags fainter) and extinction. So if the absolute 
magnitude of the source is already known, the extinction 
can be determined from the reddening and hence the dis- 
tance is directly determined (see formulae in § |3.2| ). When 
there are a number stars of the same spectral class (i.e 
approximately the same absolute magnitude and color) 
but at different distances from the Sun, these will be sit- 
uated along a line on the CMD. Increasing the distance 
shifts the star to a fainter magnitude, while extinction by 
itself shifts the stars diagonally. Hence, combining these 
two effects would mean that sources from a single spectral 
class will form a line (or stripe) on a CMD. In areas of 
the Galaxy where extinction is low this line will be almost 
vertical whereas for lines of sight into the inner Galaxy, 
the line will have a considerable angle as the extinction 
becomes the dominant effect. 

On a color - absolute magnitude diagram (CAMD), 
the majority of the stars lie in the main sequence or giant 



branch. When there is a wide range of distances being 
sampled, each part of the CAMD produces its own stripe 
on a CMD and so the overall effect is more of a broad 
smear rather than a well defined line. By far the most 
common type of giant, however, are the early K giants 
or red clump giants. These produce a very well defined 
grouping, with a density far higher than the surrounding 
areas, on a visible CAMD (e.g. Hoeg et al 1997). On an 
IR CAMD (e.g Cohen et al. 2000) the grouping is even 
tighter. For V - [8.3] CAMD the total spread of the red 
clump is about 0.5 in color and 1 magnitude in absolute 
magnitude. There are sufficient of these sources to produce 
a clearly recognizable feature in a CMD, and this was 
used in Hammersley et al. (2000) to examine lines of sight 
into the inner Galaxy. The luminosity function derived in 
Hammersley et al. (2000) shows that at the peak of the 
red clump there are about 10 times the numbers of sources 
as there are either 0.8 magnitudes brighter or fainter, and 
the sigma in absolute magnitude is about 0.3 magnitudes. 
The Mx of the center of the red clump is —1.65, although 
there is a possible error of about 0.1 magnitudes. Using 
the tables in Wainscoat et al. (1992) an Mjc — —1-65 
corresponds to a K2III and hence the J-K color is ( J — 
K) = 0.75. 

Figure [l] shows the CMD for two of the fields used 
in this study. The red clump stripe can be clearly seen. 
In order to isolate the red clump sources, we have made 
histograms of horizontal cuts (i.e. in color) through the 
CMDs, varying the vhk range used. A gaussian was then 
fitted to histogram (i.e. the region containing the red 
clump giants) to determine the position of the peak in 
each horizontal cut. Stars lying well to the left of the max- 
ima are almost totally dwarfs, as there are very few giants 
with colors bluer than the red clump stars. Stars lying to 
the right will be predominantly the M giants and AGBs, 
at least brighter than uik =13. 

In order to simplify the parameter fitting it will be as- 
sumed that the red clump giants all have Mk = —1.65 
(i.e. luminosity function becomes a Dirac delta function) 
instead of a more realistic case with some dispersion (see 
in §3.3 a discussion about this approximation). The peak 
of the histogram is used to trace the red clump giants em- 
pirically. However, it can be seen in Fig. [j] that at fainter 
magnitudes the effect of the K dwarf population becomes 
significant and swamps the red clump sources. Hence, this 
limits the method to rriK < 13.0. Figs. [jjshow the position 
of the peaks on each of the CMDs. Using the fitted trace 
for the red clump giants extinction along the line of sight 
for each field can be determined (see §|3.2j ). 

The red clump giants can be directly extracted using 
the fitted trace. Sources with a (J — K) within 0.2 mag 
of the center of the fitted red clump line were extracted, 
and the star counts were then obtained. The width was 
determined after a number of trials and is a compromise 
between the necessity to avoid other stars types but in- 
clude most of the red clump giants. If too small a width 
were chosen then any error in the determination of the 
position peak would lead to a loss of sources at that point 
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and so reduce the measured density. Using this width, 
most of the red clump giants for low extinction regions will 
be extracted. However, where the extinction is somewhat 
patchy (for instance in negative latitudes) the strip broad- 
ens and a significant number of stars are lost. Therefore, 
those regions have been excluded from the analysis. 
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Fig. 1. Color-magnitude diagrams for two studied fields. 
The higher extinction causes a slight broadening in the 
K-giants strip for in-plane regions. Maxima obtained by a 
gaussian fitting to star distribution in the red clump strip. 
The solid line shows the fitted trace that we assign to red 
clump giants population. The dashed lines show the limits 
for the red clump giants extraction within a width of 0.4 
mag. 



One of the main strengths of this method is that it 
is empirical and the only assumption being made is that 
the absolute magnitudes of all of the sources being ex- 
tracted is fixed. If, however, we had taken a theoretical 
trace for the red clump giants using a model and isolated 
the stars around this, then any difference between the pre- 
dicted and actual extinction would directly lead to errors 
in the extracted counts. Another advantage of this method 



is that the spatial information is obtained directly from 
the counts. The position of the sources can be determined 
directly from the magnitude and extinction without the 
need to resort to models. 

The error analysis presented in the following sections 
will only include the random errors in the fitting and not 
the systematic errors outlined in this section. The error 
in the determination of the absolute magnitude is taken 
to be 0.1 magnitudes, which corresponds to a distance er- 
ror of 5%. The possible error in (J — K) only affects the 
calculation of the extinction, as the peaks of the color his- 
tograms are derived directly from the data. An error of 
0.1 mag in (J — K)o would lead to an error of 0.07 mag 
in the extinction [see eq. (||)] and hence an error of 3% in 
distance. The errors in assuming that the luminosity is a 
delta function are also of the order 5%. To these should 
be added the errors in the 2MASS zero points or the as- 
sumed extinction coefficients although these are likely to 
be small. The photometric error in 2MASS sources with 
mx < 13 is ~ 0.03 mag., however the error is much lower 
when using the mean magnitude in a bin with tens of stars. 
Similarly, the error in the 2MASS-color is negligible for of 
the same reason. Hence, as well as the stated error there 
will be a systematic error in the distances, but this should 
be under 10%. 

3.2. Stars density. Disk scale-length 

Using the red clump star counts, we can directly obtain the 
star density (stars pc -3 ) by means of the stellar statistics 
equation: 



A(m) = w / r 2 D(r)4>(M)dr, 



(7) 



where A(m) is the number of stars per unit area of solid 
angle w at m in interval dm, and 4>(M) is the luminosity 
function. The bins used in this paper have dm = 0.1 in 
the range m/- = 8.5 to m/- = 13.0. 

As has been already discussed, we have assumed 
Mk=— 1-65 and an intrinsic color (J — if)o=0.75 for a 
red clump giants. The extinction Ak {nik), to a distance 
r, can be determined by tracing how (J — K) of the peak 
of the red clump counts changes with mx ■ The extinction 
is calculated for any using the measured (J-K) of the 
peak, the intrinsic mean color of the stars in the apparent 
magnitude bin (J — K )o, and the color excess definition 
and the interstellar extinction values for Aj and Ak, as 
given by Rieke & Lebofsky (1985): 

(J-K)-(J-K) 



A 



K 



1.52 



(8) 



Since we have assumed that red clump giants have 
an absolute magnitude of Mk=— 1.65, 4>(M) is replaced 
by a delta function. Furthermore, by choosing apparent- 
magnitude intervals of 0.1 mag (5m) and the related dis- 
tance intervals (Sr), the star density can be directly deter- 
mined using: 



10" 



-M k +5-A k (t) 



(9) 
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D(r) 



A(m)Srr 
wr 2 Sr 



(10) 



The galactocentric mean distance of the stars per bin 
(R) is calculated using the measured mean distance (r) in 
the following equation: 



R = J Rq + (r cos 6) 2 — 2r cos bR cos I 



(11) 



With this method, we are inverting the star counts 
of the red clump population to obtain the stellar density 
along the line of sight. This is a simple and direct method 
and although some assumptions and approximates are re- 
quired, these can be justified. Another method of inver- 
sion of star counts is presented in Lopez-Corredoira et al. 
(2000) which solved for the both density and luminosity 
function. However, this method only works for the bulge 
and would not be appropriated for the disc. 

3.3. Testing the method 

In order to to verify that the proposed method is correct 
a number of tests were made: 

3.3.1. The difference between the real red clump 
distribution and a Dirac's delta 

To further explore the form of red clump, we have cal- 
culated the distribution of K magnitudes and (J — K)q 
colors using the "SKY" model, but updating the density 
function of the giant populations (Wainscoat et al. 1992; 
Cohen M., private communication). In Fig. |], the part of 
the luminosity function pertaining to the G5 to MO giants 
is shown. The right hand plot shows the predicted form of 
the (J — K) histogram at tbk-=+11 (but with no extinc- 
tion) along a line of sight towards the anti-center. In both 
plots the predominance of the red clump giants is evident. 

In order to determine the effect of using a Dirac delta 
as the luminosity function for the red clump, we have car- 
ried out a simple calculation. The red clump giants are 
assumed to have a gaussian dispersion in magnitude with 
a = 0.3. From eqs. (0) and (0), the density in the direc- 
tion of the anticenter can then be written as: 

p'(R-R Q )= ' J 



ln(10)%/2^cr (R-Rq) 1 



dr r 2 p(r)e 



I (flog 



(12) 



where p' is the density derived with the method developed 
in this section (see § |3.2| ), while p is the real density dis- 
tribution. Using eqn. (|6j) with H = 2.0 kpc, p Q — 1, gives 
the result plotted in Fig. [| The difference between the 
two curves is the error introduced by the approximation, 
hence it is negligible. The scale-length that is determined 
assuming Mr- = —1.65 for all of the sources is H' = 1.96 
pc, i.e. the error in the calculation of the scale-length is 
around 2%. Therefore, the simplification is justified and it 
does not produce significative differences. 
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Fig. 2. Distribution of K absolute magnitudes and (J — 
K)q colors at tuk = 11 along a line of sight towards 
the anti-center corresponding to the sum of all giant 
types in the disc according to the updated "SKY" model 
(Wainscoat et al. 1992; M. Cohen, private communica- 
tion). The maxima of these plots corresponds the red 
clump, which are predominant in the giant population. 



3.3.2. Can the cut-off be masked by contamination of 
low luminosity giants with the same color than 
the red clump? 

One critical aspect of the method is the assumption that 
there are only red clump sources being extracted (i.e. 
Mk = —1-65). However, the CAMD presented in Cohen 
et al (2000, their Fig. 4), shows that there are lower lu- 
minosity sources with the same color as the red clump. 
(The range of (J — K) used for extracting the red clump 
corresponds to a V— [8.3] of about 1.2). The magnitude 
limit for the CAMD is about m\g.s\ = 6.5 and so in the 
100 pc plot of Fig. 4 in Cohen et al. (2000), the source 
count is complete for all absolute magnitudes brighter 
than 1.5. There are some higher luminosity sources with 
the same color and the red clump, but these are relatively 
few and can be ignored. However, there are significant 
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Fig. 3. Comparison between a real density p = e~~H and 
the convolved distribution of eqn. ( |l2| ) obtained under the 
assumption of Dirac delta luminosity function instead of 
the real gaussian distribution with a = 0.3. 

numbers of lower luminosity giants. These giants will have 
+2 > Mk > — 1.4 and could mask a cut off. 

In order to gauge the effect of these sources on the 
method, a simple luminosity function was built which had 
half of the giants with -1.45 > Mk > -1.85 and the 
other half of giants with +2 > Mk > —1.4. This is prob- 
ably an overestimate of the importance of the lower lumi- 
nosity sources and so should be considered a worse case 
scenario. We have carried out the same type of convolution 
in the anticenter of the previous subsubsection but with 
this luminosity function, and we have used a truncated 
exponential disc. The result is shown in |J. 

As can be observed, the only significant change is in 
the amplitude because half of the stars with Mk — —1.65 
have been lost. However, this is not important as the total 
number of the giants is not used in the analysis. Also the 
derived scale-length changes by 7%. The cut-off, however, 
is not masked as there is a clear change of slope in logp 
around d w 5 kpc. Hence, if there were a cut-off at about 
13 kpc it should still be visible even with the presence of 
possible lower luminosity giants. As will it be shown below 
(Figs. |7[ ^), we do not observe any change of slope in log p, 
so the absence of cut-off seems clear and cannot be criti- 
cized as an effect of the contamination of low luminosity 
giants. 

3.3.3. Dwarf contamination 

Another critical aspect that must be considered is the 
dwarf contamination. There is an important contribution 
to the stars counts from dwarf stars for faint apparent 
magnitudes (vhk > 14; Figs. |l|). The selection method 



Fig. 4. Comparison between a real density p = e~"R , trun- 
cated at R = 12.9 kpc (distance from the Sun: 5 kpc) 
and and the convolved distribution obtained under the as- 
sumption of Dirac delta luminosity function instead of the 
real distribution with 50% of the stars with Mk = —1.65 
and 50% of the stars between +2 > Mk > —1.4 (lower 
luminosity giants). 

used here cannot separate between these two star popula- 
tions at these magnitudes, as the main sequence and the 
K-giant strip overlap. In order to estimate the dwarf con- 
tamination for rriK < 13.0, the luminosity and densities 
functions of the giants and dwarf stars from the model 
of Wainscoat et al. (1992) were used to calculate the pre- 
dicted ratio between dwarf and the total of giants plus 
dwarfs in the counts. The result is shown in Fig. |^. It 
is clear that for faint apparent magnitudes [jtik — 14) 
the dwarf contribution is important, accounting for al- 
most 50% of the total counts. However, the contribution 
decreases at brighter magnitudes such that for m,K < 13 
less than ~ 10% of the stars are dwarfs. This was cal- 
culated for off-plane regions, where the extinction is al- 
most negligible. It is expected that for regions close to the 
Galactic plane, where the extinction is higher, the dwarf 
contamination will be even less, because the reddening will 
shift the giants to the right whilst the dwarfs will be less 
affected, and so the dwarf-giant separation is improved. 
Hence, here we will only use the red clump star counts for 
ttik < 13.0, where the dwarf contamination is small. 

3.3.4. The ambiguity between reddening and 
metallicity effects for the clump giants 

When the metallicity decreases the stars become bluer, 
and hence the reddening could be underestimated. 

We have assumed that the stellar populations through- 
out the disc is homogeneous and therefore, there are 
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Fig. 5. Percentage of predicted dwarf stars in the total 
star counts extracted from the color-magnitude diagram 
for I = 180°, b = 0°. 



no large gradients in the metallicity. However, even in 
the case of small metallicity variations (Cameron 1985, 
has measured a radial metallicity gradient of -0.11 dex 
kpc -1 ), the change to the result is small as the red clump 
color is practically independent of age and metallicity 
(Sarajedini et al. 1995). The variations are limited within 
A(B — V)o < 0.05 (Sarajedini et al. 1995), which means 
A( J— K)q < 0.022, i.e. an error in the extinction less than 
0.015, which is negligible. 



tocentric distance before reaching the outer disc, which 
would increase the error in obtaining the scale-length. 

The main characteristics of the chosen fields are sum- 
marized in Table [l]. Column (1) and (2) give the Galactic 
coordinates of the center of the field. Column (3) is the 
total area in the field which has 2MASS data. The total 
number of stars with < 14.5 is given in column (4), 
and finally column (5) gives the total number of red clump 
giant stars extracted from the color-magnitude diagram, 
with nik < 13. 

3.5. Results 

A sequential fitting procedure has been used to obtain 
the model parameters for eqs. ([!]) and (||) (or eqs. (||) and 
(^)). First, to restrict the number of free parameters we 
have used only the data at b = to obtain the radial 
scale- length H (i.e., eq. (^|)), so that there would be no 
dependence on height. This scale-length is then used in 
the off-plane regions to derive the values of the scale-height 
(h z (R)) and flare (h R . flare ). 

An exponential disc (eqn. |6|) was fitted to the data 
with rriK < 13.0 in the plane for longitudes of 155°, 165°, 
180° and 220°. The value of H which gives the minimum 
X 2 (x 2 = 143.2 for a number of data points N = 184, see 
Fig. p|) and assuming that the errors for each data point 
are poissonian was for: 



H = kpc. 



(13) 



300 



3.4. Selected regions 

Twelve fields were chosen from the available 2MASS data, 
which are summarized in Table [l| In each field a large area 
of sky was covered (between 2 and 5 square degrees) as a 
large number of stars are needed to determined the color of 
the peak of the red clump distribution. The high number 
of red clump giants also minimizes the poissonian noise in 
the star counts and so improves the subsequent fitting of 
the density distributions. The fields were 5° in longitude 
by 1° in latitude, although the actual area covered depends 
on the available 2MASS data within each region. The star 
density is changing far more rapidly in latitude than lon- 
gitude so it is important not to include a wide spread 
of latitudes in each region. Positive latitude regions were 
selected because these are less affected by patchy extinc- 
tion, which causes a significant dispersion of the red clump 
strips. We have only used regions 150° < |Z| < 225° so as 
to avoid the northern and southern warps, which produce 
large asymmetries between positive and negative latitudes 
(see §||). Furthermore, the lines of sight in 45° < \l\ < 90° 
travel for large distances at approximately the same galac- 
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Fig. 6. x 2 vs. H of the fit of red clump star density in the 
plane of the model dq) . 



As described in 53.1 



the red clump sources are ex- 
tracted from the general star counts by means of the 
(J — K) CMD. Inserting these number counts in eq. (||), 
with the help of eq. (||) with the calculated value for H, we 
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Table 1. Selected regions used with the red clump extraction. 
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derive a density for these sources in the solar neighborhood 
p©,red clump = 1-31 x 10 -5 stars pc~ 3 . The disc model with 
the optimum value of H is shown in Figure Q together with 
the density distribution of the red clump stars derived as 
explained in §3.2. The small differences in the density for 
the different lines of sight stem presumably from the dif- 
ferences in the patchiness of the extinction, which spreads 
out the red clump giants, or contamination. Even so, it 
is noticeable how well the analytical function matches the 
data. 

Using the value of H determined above with off-plane 
regions at I — 180° and I = 220°, allows the disc scale- 
height and the flare scale-length to be examined. As be- 
fore, the disc model of eqns. (||) and (0) are matched to 
the number counts by minimising the x 2 of the fit, being 
now the free parameters h z {R®) and hRji are . The best fit 
(X 2 = 245.4 for N = 460) is for: 

h z (R Q ) = 3101S PC, 



The errors are calculated for H = 2.10 kpc and other 
values of H could slightly change the error ranges. The 
stated errors are determined to be changing the param- 
eters such that x 2 becomes a factor / larger than the 
minimum (which corresponds to a 68% confidence level 
and is specified in statistics manuals). The correlations 
between the different parameters do not significatively 
change the errors. Figures || and ^| show the comparison of 
this model with the observational data at different heights 
above the plane. Since the heliocentric distance to the ex- 
tracted stars can be measured (see § |3.2| ), their distance 
from the Galactic plane is also known. Both figures show 
that there is an increase in scale-height with increasing R. 
Furthermore, in Fig. ||, the data show a slower decrease of 
stellar density with z at larger R (h z = 580 pc, 1055 pc, 
1905 pc for R = 10, 12, 14 kpc respectively). Both of these 
effects can only be produced by a flare. 

This result cannot be produced by a combination of 
a thick disc and a thin disc with constant scale-heights 
and the scale-lengths because, then the global scale-height 
would be independent of R. If there is a thick disc with, 
for example, h z = 910 pc, a density on the plane of 5.9% of 
the thin disc (Buser et al. 1999) and assuming that there 
is not significative difference in the abundances in the red 
clump population, then the thick disc would become dom- 
inant for z > 1300 pc whereas for z < 500 pc the thin 
disc would be by far the more important. Fig. || shows 
very different slopes for different R at z < 500 pc, and 
even the combination of thin and thick disc with differ- 
ent scale-lengths cannot reproduce this effect because the 
thin disc is highly dominant for z < 500 pc. Furthermore, 
the thick disc cannot produce significant errors in the cal- 
culated scale-length, since this was derived in the plane 
where the thin disc always dominates. The scale-height, 
however, would be representative of an average of thin 
and thick discs. The average source height in the data set 
used is ~ 300 pc, where the thin disc is ~ 10 times denser 
than the Buser thick disc. 

Another noticeable feature of Figures [?] and ^ is that 
the density follows the exponential to at least 15 kpc. As 
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has already been discussed, if the disc had a sharp cut- 
off then this would cause an abrupt drop in the density. 
Hence the absence of such a sharp drop implies that there 
is no abrupt cut-off to the disc until at least ~ 15 kpc. By 
15 kpc the star density has become too low to determine 
if there is a cut-off beyond that point. 
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Fig. 9. The red clump star density obtained from the data 
and the lines showing a model (Q) with H = 2.1 kpc, 
h z (R(7))—310 pc, hnjiare = 3.4 kpc at I = 220° and the 
latitudes stated on the plot. 



Fig. 8. The red clump star density obtained from the data 
and the lines showing a model (||) with H = 2.1 kpc, 
/i z (i? Q )=310 pc, hjiji are = 3.4 kpc. Note that the data 
show a slower decrease of density for larger R, and this fact 
can be modeled with the flare h z = 580 pc, 1055 pc, 1905 
pc respectively for R = 10, 12, 14 kpc (the data belong to 
a range of 0.5 kpc centered in these galactocentric radii). 



3.6. Flare in the inner Galaxy? 

The regions used in the previous subsection all had |Z| > 
90°, and hence the lines of sight run external to the solar 
circle (R > R@). Therefore, the results are valid only for 
the outer disc. Moreover, the flare scale-length of around 
3.4 kpc is really only a description of the flare for R > 10 
kpc, as between 8 and 10 kpc there is little data which is 
also very noisy (Fig. [?]). Hence, whether this law is appli- 
cable to R < 10 kpc is something that cannot be deduced 
from the previous analysis, and data with a range of lower 
galactocentric distances is required. 

In order to examine the effect of the flare in the in- 
ner Galaxy this method should not be applied to lines 
of sight which cover only small a range of Galactocentric 
distances inside the solar circle, as these provide only a 
small range of values with which to attempt fitting R. 
For instance, on plane regions with 45° < I < 90° would 
only give values of R between 6 and 8 kpc. Therefore, 
in order to explore the inner disc using this method, we 
must use regions with |Z| < 45°. Furthermore, the near 
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Fig. 10. The red clump star density obtained from the 
data and the lines showing a model (S) with H — 2.1 kpc, 
Ii z (Rq)=310 pc and different Hare scale-lengths (fiRjiare) 
at / = 30°, b — 6°. As can be seen, a flare is necessary 
in the inner Galaxy, and its best scale-length is around 9 
kpc. 

plane regions with |/| < 30° should also be avoided, as 
the inner Galaxy components become important, e.g. the 
in-plane bar (Hammersley et al. 1994, Lopez-Corredoira 
et al. 2001). At present the available 2MASS data in the 
regions of interest is very restricted and so we will concen- 
trate in only one region: I = 30°, b — 6°. 
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Figure [H] shows the density of the red clump giants 
along this line of sight determined using the same method 
as in §3.2. As well as the measured density, various model 



predictions are plotted. It shows that a flare is necessary to 
describe the data with a disc of H = 2.1 kpc and h z (R Q ) = 
310 pc. However, the scale-length of the flare cannot be 
as low as 3.4 kpc, but rather around 9 kpc. The scale- 
height in the inner Galaxy is smaller than the scale-height 
in the solar neighborhood, however the difference is only 
h z (5 kpc) w 0.73h z (R Q ). A flare scale-length of 3.4 kpc 
would imply h z (5 kpc) w 0.43/i z (i2 Q )). Hence, although 
the flare begins well inside the solar circle, the increase of 
the scale-height with R is significantly slower in the inner 
Galaxy. 

This result is consistent with previous work. For in- 
stance, Kent et al. (1991) has measured h z (5 kpc) = 
O.67/i z (i?0) by fitting the parameters of the disc in a 2.4 
/jm map of the northern Galactic plane. 



Eaton et al. (1984) luminosity function 




4. Using the star counts to obtain the parameters 
of the disc 

The second method used for determining the parameters 
for the stellar density is by the fitting of star counts in 
different regions. For each of the line of sight centered 
on Galactic coordinates (l,b)i, (i is the field number) the 
cumulative star counts (Nk) observed in K, up to a mag- 
nitude rriK in a given solid angle uj is the sum of the stars 
within the beam brighter than the given apparent magni- 
tude (Bahcall 1986). Assuming that the luminosity func- 
tion does not vary with the spatial position, this is 

N K (iriK,l,b) = lo 

x / <S>K(M K (r,l,b))p(R(r,l,b),z(r,l,b)ydr, (15) 
Jo 

where 

M K (r,l,b) =m K + 5 - 5 log 10 r - A K (r,l,b), (16) 

i-Mk 

$k(M k )= / <j> K [M)dM, (17) 

J — oo 

<Pk is the normalized luminosity function for the K band 
in the disc and Ax(r, I, b) is the K extinction along the 
line of sight (I, b) up to distance r. We have used a stan- 
dard disc luminosity function given in Eaton et al. (1984, 
see Fig. |ll|) and the extinction model given by Wainscoat 
et al. (1992) (except with hn,dust = 3.2 kpc to scale the 
distance to the center of the Galaxy to 7.9 kpc). Fitting 
could also be used to constrain the extinction, however it 
is preferable to have as few free parameters as possible. 
The extinction model of Wainscoat et al. (1992) is ax- 
isymmetric and does not include a flare or a warp in the 
dust distribution. However, K-b&nd counts are not very 
sensitive to small variations of the extinction and were ex- 
tinction the goal of the paper shorter wavelengths would 
be more appropriate. Therefore, as the lines of sight are 



Fig. 11. K-band normalized luminosity function used in 
the second method of this paper, from Eaton et al. (1984). 



in general looking towards the outer Galaxy and the ex- 
tinction is low, a simple extinction model will suffice. 

The stellar density (p) parameters (eqn. ||) were 
determined using the available 2MASS data at b = 
0,±3,±6,±9 degrees; 45° < I < 315°. On or near-plane 
regions were preferred because these are more sensitive 
to the variations in the parameters, and there are signifi- 
cantly more stars. Again the central regions of the Galaxy 
were avoided so that other components other than the disc 
would not be included. It should be noted that even when 
the local stars dominate the counts there are still sufficient 
giants to constrain the large scale structure. Each region 
used an area of 1 square degree (Al = 1°, Ab = 1°). 
However, some areas are only partially covered and in 
these cases the counts were divided by the actual area 
available to obtain the counts per unit area. In total 820 
lines of sight with areas between 0.5 and 1.0 deg 2 are used. 
This gives a very large area coverage along the Galactic 
plane, and so can constrain the parameters of the disc. 
Other authors have used 2MASS star counts to study the 
disc however the number of regions used is much lower. 
Ojha (2001) used seven independent regions quite far from 
the plane but does not study the warp, flare and cut-off 
whilst Alard (2000) only looked at three strips in latitude 
across the plane. 

As noted, this method assumes that all sources have 
the same spatial distribution, i.e. the scale-height and 
scale-lengths in a given galactocentric radius are the same. 
For the scale-height, in particular, this is known not to 
be true with the younger components having a scale- 
height far less than the older components. However, the 
young components only make a significant contribution 
to the counts at the brightest magnitudes (rn^c < 8) and 
even then only on the plane towards the inner Galaxy 
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(Hammersley et al. 1999). For this study we have used 
2MASS counts up to tvik — 14.0. This magnitude limit, 
coupled with using lines of sight away from the inner 
Galaxy, means that the counts will be completely dom- 
inated by the old population of the disc rather than the 
young disc and spiral arms. Hence the assumption that all 
of the sources will have, at least, very similar scale-heights 
is justified. 

We have limited the maximum apparent magnitude in 
K band to 14.0. This is well below the limiting magnitude 
of 2MASS and ensures that there are complete statistics in 
all regions. Moreover, with this limiting magnitude, stars 
brighter than Mk — can be detected to over 6 kpc 
which means the results are sensitive to the shape of the 
very outer disc. However, one of the curious features of 
stars counts is that the fainter the magnitude the closer 
on average are the sources being detected. This is caused 
by the steeply rising luminosity function and means that 
the dwarfs will dominate the fainter counts. The result is 
that the red clump giant method, which isolates a specific 
group of stars, will provide a better contrast for the more 
external part of the disc. 

As discussed previously, the "old disc" analyzed here 
includes all the old stars of the disc within the selected re- 
gions (average latitude 4-5 degrees). If there is a thick disc 
like the one proposed by Buser et al. (1999), the mean con- 
tribution to the star counts in the selected regions would 
be around 10%, assuming a thick disc luminosity func- 
tion similar to that of the halo in Wainscoat et al. (1992). 
Hence, the scale-height determined here would be repre- 
sentative of an average of the thin and thick discs but with 
90% of the sources coming from thin disc. 

Star counts cannot easily distinguish between the cut- 
off and the loss of stars in the plane due to the flare. 
However, it was shown in the previous section that it there 
is no cut-off to at least 15 kpc, and so we will continue to 
assume there is not a cut-off. 

4.1. Results 

As in §||, the fitting of the parameters will be in two steps. 
The first step is to fit H using the on plane regions (b = 0) 
and then fit h z (Ro) and /ifl,flare using all of the regions. 

Using the available regions in the plane the minimum 
X 2 ( x 2 = 452, N = 114) was obtained for (see Fig. |l|): 

H = 1.9li°;*> kpc. (18) 

The stellar density in the solar neighborhood is p© = 
5.3 x 10~ 2 star pc~ 3 . This is somewhat lower than the 
value measured in other surveys (e.g., 0.13 star/pc 3 , 
Bahcall 1986). This discrepancy is most probably due to 
using different luminosity functions, here we used the nor- 
malized luminosity function by Eaton et al. (1984). In any 
case the normalization does not affect the other disc pa- 
rameters. The best fit is shown in the figure [l3| (dashed 
line) . 
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Fig. 12. x 2 of the fit of star counts in the plane for the 
model (||). 
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Fig. 13. Fit of star counts in the plane. The dashed line 
is the best non-warped model (H — 1.91 kpc). Solid line 
is the best warped model (H = 1.97 kpc; e w — 5.25, 4> w — 
-5°, A w = 2.1 x 10~ 19 pc). 

The other two parameters were then fitted using 
all the available data at 45° < / < 315° and b — 
o ,±3°,±6 o ,±9°. The minimum X 2 is found for 

h z (R Q ) = 300t\l pc, 

hit ,flare — 

4.6 ±0.5 kpc. (19) 

These errors are for H = 1.91 kpc and other values 
of H could slightly change the error ranges. In this case, 
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the best fit was for x 2 = 3233 with N = 820. The fits are 
shown in Figures [l^ and [l4| (dashed line). 

In general there is a good fit, although there are some 
asymmetries between 45° < I < 180° and 180° < I < 315° 
and also between b > and b < 0. These can be attributed 
principally to the warp (see §|^), but there are also some 
other minor fluctuations, possibly due to the patchiness 
of the extinction or contamination from other Galactic 
components. 

5. In search of the stellar warp 

If there were neither a warp nor other asymmetries, the 
ratio of counts above and below the plane (Fig. |l5|) should 
be nearly the same for all Galactic longitudes, with only 
a small offset due to the height of the Sun above the 
plane. However, the plot clearly show a sinusoidal behav- 
ior. There is an excess of counts at positive latitudes with 
respect the negatives for I < 180° and the opposite for 
I > 180°. A warp, together with a small offset due the 
height of the Sun above the plane, can explain what is 
seen. 

A model similar to eqn. ([|) will be used, but with 
| z — z w \ instead of \z\, where z w is the average elevation 
of the disc above a plane which is parallel with b = but 
runs though the local center of mass (i.e. 15 pc below the 
Sun). 

z w = [C w R{pcf w sin(</> - <f> w ) + 15] pc. (20) 

Hence, the model describes the warp as a series of tilted 
rings, where the amount of tilt is the galactocentric dis- 
tance raised to a power e w . The 15 pc term is due to 
the height of the Sun above the plane (Hammersley et al. 
1995). Since the value of about 15 pc is well constrained 
by other studies we will not carry out an analysis of this 
parameter. The two warp parameters, e w , and the galac- 
tocentric angle 4> w {4> w q = 0), are determined from the 
division of star counts data in Fig. [l^. Once the parame- 
ters are fixed then C w is obtained from: 

C w = ) (21) 

where di = ( — lj with N(m,b) the observational 
cumulative counts, cr; the errors of di, and ij is the model 
prediction for di with C w = 1. 

In this way, there are only two free parameters in 
the model and so it is simpler to constrain the warp. 
The model assumes that the northern and southern warp 
are symmetric, apart from the effect of the height of the 
Sun above the plane. This is an approximation and it is 
known that the gas warp is not symmetric (Burton 1988). 
However, the asymmetries of the warp are produced over 
galactocentric distances R « 13 kpc (distance from the 
Sun: r > 8 kpc for I — 250°), where there are few stars 
from the total sample of stars with < 14.0. 
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Fig. 15. Positive counts divided by the negative counts. 
The solid line stands for a warped model with H = 1.91 
kpc, h z (R e ) = 300 pc, tiRjiare = 4.60 kpc; e w = 5.25, 
<f>w — —5°, A w = 2.1 x 10 -19 pc. The dashed line stands for 
the best "tilted disc" model with H = 1.91 kpc, h z (R Q ) = 
300 pc, tiRji are = 4.60 kpc; 6 C orroctcd = b + 1.64° sin(Z + 
170°). 
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5.1. Results 

The best fit is found for: 
e w = 5.25 ±0.5 

4>w — — 5 ± 5 degrees, (22) 

which gives C w = 2.1 x 1CT 19 pc with X 2 = 332.7, 
N = 217. The best fit is shown in Fig. |l5| (solid line). 
Although there are few points which do not fit very well 
(perhaps due to irregularities in the extinction or in the 
distribution, contamination of the spiral arms, etc.) again 
the general fit is good. Hence, a warp does accurately de- 
scribe the deviations from the best symmetric disc. In Fig. 
|lq, a representation of the warp is plotted. 
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Fig. 16. A schematic of the warp showing the relative lo- 
cation of the Sun. The rings arc each kpc, and the line of 
the maximum warp is shown. 



5.2. Could this be a local tilted disc effect? 

An alternative hypothesis that could explain the asym- 
metries is that the disc is locally tilted with respect 
to the global galactic plane. Hammersley et al. (1995) 
showed that a similar asymmetry COBE-DIRBE surface 
brightness-maps could be caused by a local tilt to the disc. 
However, Hammersley et al. (1995) deliberately chose re- 
gions well away from the plane, so that the results would 
be dominated by local sources and there would be little 
contribution from regions where the warp becomes sig- 
nificant (R >11 kpc). Star count data provide direct in- 
formation on the distance to the sources, something not 
available with surface brightness-maps and so can better 
distinguish between a tilt and warp. 

In Fig. [l5], it is shown the best fit for a tilted disc 
model, i.e. assuming that the asymmetries are only due to 
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the height of the Sun above the plane and the local tilt 
of the disc plane, which corrects the Galactic latitude as 
(Hammersley et al. 1995): 

z w = 15 pc, 

^corrected = b + a t sm(l + l t ). (23) 

The best fit for the tilt corresponds to 

It = 170 ± 20 deg, (24) 

which gives a value of at = 1.64° and % 2 = 517 (N — 
217). The orientation found is in good agreement with 
Hammersley et al. (1995), but the tilt is some 4 times 
greater indicating sources at a significantly larger distance 
are being detected. 

The predictions for the tilt are significantly poorer 
than for the model of the warp, as can be seen in Fig. 
|l5| (dashed line). Hence, whilst a small tilt is sufficient 
to describe the geometry of local disc, beyond a few kpc 
this simplification is no longer valid. Interestingly the av- 
erage tilt predicted by the warp model over the distance 
range used in Hammersley et al. (1995) (i.e. about 2 kpc) 
is still about twice as large as that measured from the sur- 
face brightness map. This suggests that the warp actually 
starts just outside the solar circle, although it should be 
noted that the predicted effect in the solar vicinity is small 
in any case and the star counts are not sensitive for very 
short distances from the Sun. 

5.3. Confirmation of the warp with the red clump 
giants 

A warp should also leave a clear imprint in the distri- 
bution of the red clump giants (^). Since the warp will 
affect the distant sources most, the asymmetries should be 
greater at larger galactocentric distances. This can be seen 
in Fig. [l7|, where the ratio of densities at positive and neg- 
ative latitude is shown for the region I = 220°, b = ±6°. 
The data show a strong asymmetry with most points hav- 
ing a ratio < 0.8 whilst for those beyond 13 kpc it drops 
to under 0.6. This again shows that the asymmetries in 
the counts are larger for distant sources. To confirm that 
asymmetries cannot be due to the height of the Sun above 
the plane, the dashed line in Fig. [ll] shows the predicted 
level of this effect. The predicted warp, however, is shown 
as the solid line and this does provide a reasonable fit. 

Another aspect which needs to be commented upon is 
that the measured flare scale-length may be affected by 
the warp, even in regions where the warp is not very con- 
spicuous (I — 220°). Fig. p7| shows that the ratio may be 
as low as 0.6, i.e. an error of around 20% for each of the 
densities at positive and negative latitude. This is one rea- 
son why the results from (jlj) should be taken with some 
care. However, in the second method, the relative contri- 
bution from the very distant stars (R >^ 12 kpc) is much 
lower and the effect of the warp is smaller (compare the 
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Fig. 17. Ratio of densities between positive and negative 
latitude for (/ = 220°, b = ±6°). The solid line shows 
the predicted ratio for the warp model with H = 1.91 
kpc, h z (R Q ) = 300 pc, hujiare = 4.60 kpc; e w = 5.25, 
(f) w = — 5°, C w = 2.1 x 10~ 19 pc. The dashed line shows 
the effect of the height of the Sun above the plane (+15 
pc). 



results between §4.1 and §p.5|). Hence, the parameters de- 



rived from the star counts are probably more trustworthy 
than those derived from the red clump giants in §g. 

5.4. Coincidence of stellar and gas warp 

The measurements of stellar warp presented here are to- 
tally independent of those used to determine gas warp 
(Burton 1988, Diplas & Savage 1991). It is noticeable how- 
ever, that the stellar warp is nearly identical to that of the 
gas. Not only is the angle of orientation of the maximum 
amplitude ( 4> w 85°) coincident, but the amplitude as well 
(see Fig. |lj) . The northern warp in both the gas and stars 
are very similar to R > 13kpc, as shown in Fig. Q The 
southern gas warp differs from the northern gas warp, be- 
coming constant or decreasing with height for R > 13kpc. 
However, since a large majority of the stars along the line 
of sight of southern warp (I ss 250°) have R < 13 kpc 
(r < 8 kpc), the data are not sensitive enough to detect 
the differences between the northern and southern warp. 

Other authors have also looked for the stellar warp. 
Djorgovski & Sosin (1989) obtained a similar result with 
IRAS sources. Alard (2000) also measured the stellar warp 
using 2MASS data, and also concluded that the amplitude 
is roughly coincident with the gas warp. Although the 
present results are more accurate, Alard (2000) showed 
that the amplitude of the stellar warp is consistent with 
the HI to within about 30%. Drimmel & Spergel (2001) 
found from COBE-DIRBE data that the stellar warp is 



Fig. 18. Maximum amplitude of the stellar warp (solid 
line: z w = 2.1 x 10~ 19 i?(pc) 525 pc, which is the best power- 
law fit to the data) in comparison with the one measured 
by Burton (1988) for the northern and southern warp gas 
scaled to Rq = 7.9 kpc (dashed and dot-dashed lines). 

more or less coincident with the stellar warp found here 
(and Burton(1988) gas warp): z w 25 pc/kpc 2 (i? — 
7 kpc) 2 and Freudenreich (1998) obtained similar results 
in his independent determination of warp parameters. 

From the presented results, it is clear that the form of 
the stellar warp is identical to that of the gas warp. Hence, 
whatever is causing of the warp affects both to the stars 
(including the old population) and gas in approximately 
the same way. Whilst the stellar and gas warps are coin- 
cident for the Milky way, this may not necessarily be the 
case in other galaxies. 

5.5. A new fit of the disc including the warp. 

After obtaining the parameters of the warp, which is a sec- 
ond order effect on the total counts, it is possible update 
the disc model to include the warp and then redetermine 
the parameters using the cumulative star counts. As ex- 
pected, there are no major changes with respect to the 
previous values (fill), (O): 



H = 1.97 



-0.15 
-0.12 



kpc 



hRJlare = 5.0 ± 0.5 kpC 

h z (R & ) = 285± 8 12 pc. 



(25) 



The errors in h,Rji are and Ii z {Rq) are for the fixed value 
of H — 1.97 kpc. Other values of H could change very 
slightly the error ranges. This best fit gives po = 5.5 x 10 -2 
and corresponds to \ 2 = 352 for the data in the plane 
(N = 114), and x 2 = 2049 for the whole sample (N = 
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820). The fits are plotted in Fig. |T| and |T| (solid line). 
Using eqn (^), this provides a constraint on the intrinsic 
scale-length of the disc: 



3.3±°i kpc 



(26) 



It should be remembered that these parameters have 
between determined between about R — 6 and R = 15 
kpc (between O.75i?0 and 1.9R®). No data comes from 
R < 6 kpc as the limit in longitude is 45° and beyond 
15 kpc the data becomes uncertain. 

6. Comparison of both methods for the flare and 
discussion 

In §||, red clump giants gives values compatible with those 
from the star counts except for the hnji are , which is 2.8er 
lower (3.4 kpc instead of 5.0 kpc). It is unlikely that this 
is a separate population with different flare, because the 
scale-height at the solar circle is coincident. We suspect 
the difference stems from a selection effect. Apart from 
the fact the measurement of red clump giants is affected 
by the warp a nd hence the value of 3.4 is somewhat less 
trustworthy (§ |5.3| ), it must be remembered that most of 
the giants are at R > 12 kpc, i.e. are dominated by distant 
stars. However, the flare sca le-length is not constant, it is 
higher for smaller R ( §3J3 shows that h^jiare ~ 9 kpc 
between R = 4 kpc and R = 7 kpc). Hence, a systematic 
difference between the two values should be expected. 
The flare scale-length can be expressed as: 

h-RJlare W 12 - 0.6i? kpC, 



(27) 



where R is in kpc. This gives the values of 9.0 kpc for 
R = 5 kpc (close to the average galactocentric distance of 
the sources along I — 30°, b = 6°); 5.0 kpc for R = 11.7 
kpc (representative of the range between 6 and 15 kpc); 
and 3.4 kpc for R = 14.3 kpc (representative of the range 
occupied by the distant giants). In the solar neighborhood, 
the value of hnjiare would be 7.3 kpc. 

In comparison, Alard (2000) has described the flare in 
2MASS data by h z = h z {R Q ){\ + (0.32 kpc _1 )(i?- R Q )). 
This is similar to our results in the range 8 < R < 12 kpc. 
Hence these results are in agreement with Alard (2000) 
and the thickening of the disc explains the rapid drop of 
the density in the outer Galaxy without the need to invoke 
a cut-off for R < 15 kpc. Indeed, as we have seen in §||, 
the cut-off at R < 15 kpc is not only unnecessary but 
inconsistent with the data. 

The flare, like the warp, is another characteristic of the 
gas in the outer Galaxy. Again the parameters determined 
for the stars are very similar to that of the gas (Wouterloot 
et al. 1990), although the gas scale-height is always far 
lower. 

6.1. Application of the disc model to brighter apparent 
magnitudes 

An aspect which must be examined is the possible varia- 
tion of the measured local stellar density p© when using a 



range of magnitudes other than uik < 14.0. In principle, 
the given density distribution should be independent of 
the magnitude range used, however, there may be errors 
in Eaton's normalized luminosity function that has been 
used here. Alternatively there could be source confusion at 
the faintest 2MASS magnitudes, although such an effect 
would at most only affect the sources within a magnitude 
of the survey limit in the regions measured. Both errors 
would affect the measured density, p® . Figure |l9| shows the 
variation of p® with limiting magnitude mx < rriK.max- 
This can be described as 



p e w 0.2402 - 2.642 x 10^%^,, 



+ 9.446 x 10- 4 rr4 ma:c , 



(28) 



The data were obtained by normalizing the star counts 
with a circular region centered in (I — 177°, 6 = 3°) with 
radius=2 degrees. The form of variation is not consistent 
with there being confusion and so it has to be due to errors 
in the form of Eaton's luminosity function. 

An alternative possibility for this variation could have 
been the density parameters are not correct and for this 
reason p® is not constant with m^. This can be checked 
by fi tting of the star counts with the model described in 



5.5 but changing the limiting magnitude and hence pQ 
It is found that the model still predicts 



using eqn. |2S 

the measured star counts (some excess in the anticenter 
could be due to some contaminant of bright sources) . As 
an example Fig. ^o] shows the fit for the whole plane for 
m K < 9.0. 



1=177 , b=3 



0.15 




0.05 - 



Fig. 19. Variation of the amplitude p® as a function of the 
maximum if-band magnitude when we use the disc model 
of this paper together with Eaton et al. (1984) luminosity 
function. Solid line is eq. (E9). 
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Fig. 20. Data and fit of star counts for tuk < 9.0. The 
dashed line is the best warped model (H = 1.97 kpc, 
h z (R & ) = 285 pc, h R ji are = 5.00 kpc; e w = 5.25, 4> w = 
-5°, A w = 2.1 x 10~ 19 pc). 

6.2. Application to high latitudes 

Fig. [H] shows how well the model fits the counts at high 
latitude regions. Even near polar regions (6 = 90°) the 
model predictions do not systematically deviate from the 
counts. It should be noted that the model does not include 
a bimodal distribution, containing both a thin and thick 
disc with different populations and scale- heights. None the 
less using the thin disc alone the model gives an accu- 
rate prediction of the counts, for the magnitude range be- 
ing used. It is beyond the scope of this paper to discuss 
whether this result should be interpreted as an absence of 
thick disc or whether it is simply that the analysis pre- 
sented here is not sensitive to the thick disc. A more de- 
tailed examination of high galactic latitudes, including an 
analysis of the colors and luminosity function, would be 
required. 

7. Conclusions 

In this paper we examine both quantitatively and qualita- 
tively the large-scale structures in the Galactic disc. Two 
independent methods are applied to the 2MASS data. The 
first isolates red clump giants in a number of regions near 
the plane, the counts are then inverted to directly ob- 
tain the disc density. The second fits the parameters of a 
disc model in 820 independent near-plane regions. Results 
from both methods are consistent with each other imply- 
ing that both methods are valid. For instance, magnitude 
and color of the red clump giants must be close to the cor- 
rect values otherwise there would be significant difference 
between the two methods for H and /i. z (i? Q ). 
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Fig. 21. Star counts in high latitude regions. Note how 
the model derived in the present paper (fitting low lat- 
itude regions) also fit approximately these high latitude 
regions. Data come from the available covered regions in 
the different coordinates within a radius of 2 degrees. 

The distribution of sources within the disc is well de- 
scribed by an exponential distribution in both the galac- 
tocentric distance and height above the plane. There is 
a strong flare, i.e. an increase of scale-height towards the 
outer Galaxy. This flare starts well inside the solar cir- 
cle, hence there is a decrease of the scale-height towards 
the inner Galaxy. This is consistent with the observation 
that the large majority of galaxies have flares (de Grijs & 
Peletier 1997) and theoretically consistent with some sce- 
narios of self-graviting stellar discs (Narajan & Jog 2002). 
The flare removes the need to invoke a cut-off in the outer 
Galaxy, and certainly there is no evidence for an abrupt 
cut-off in the stellar disc, to at least R < 15 kpc (=4.5/i#). 
For the analysis it was asssumed that the red clump lumi- 
nosity function could be described as a delta function and 
whilst it is shown that this could have some minor effect 
on the result, it could not mask a cut off. 

In other spiral galaxies, the question of the cut-off is 
still not very clear. Pohlcn et al. (2002) find that typical 
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radial profiles are better described by a two-slope expo- 
nential profile, characterized by an inner and outer scale- 
length separated at a break radius, rather than a sharply- 
truncated exponential model. Florido et al. (2001), how- 
ever, argue that there are some cases of abrupt truncation 
in the near infrared emission of the stellar profile. Kregel 
et al. (2002) finds that 60% of the spiral galaxies are trun- 
cated, although he also finds that the truncation is beyond 
4.4/1^ for galaxies with hn < 4 kpc; so probably, if there 
existed a cutoff in our Galaxy, it would not be detected in 
R < 4.5/ir. 

Another noticeable feature is the warp in the old stel- 
lar population whose amplitude is coincident with the 
amplitude of the gas warp. Figs. ^2] and ^3] show the 
form of the disc. Sanchez-Saavedra et al. (1990, 2002) 
or Reshetnikov & Combes (1998) claim that most spiral 
galaxies are warped at optical wavelengths, but it is not 
easy to disentangle the stellar and gaseous warp. 





-15 -10 -5 5 10 15 

R(Kpc) 



Fig. 22. Cut in the plane YZ of the Galaxy as described 
by the final disc model. 

Various parameters of the old stellar disc have been 
measured based on an analysis over galactocentric dis- 
tances between 6 kpc and 15 kpc. The study of the inner 
disc would be somewhat more complicate, especially in 
the plane, because the disc component is no longer isolated 
and it is mixed with other possible components: bulge, bar, 
ring, etc. As is shown in Lopez-Corredoira et al. (2001), 
the disc might be a Freeman type II disc and inside about 
3.5 kpc the disc would be truncated; this deserves another 
detailed study which is beyond the scope of this paper. At 
present we can talk about the distribution at R > 5—6 kpc 
and representative of the mean old disc at low latitudes 
(with mean \z\ around 300 pc). If there exists a thick disc, 
the results for higher latitudes would be substantially dif- 
ferent. In the following are presented the results from the 
star counts, since they are potentially more accurate. 



Fig. 23. Isodensity surface with p — 0.005 star/pc 3 as 
described by the final disc model. Vertical scale multiplied 
by a factor 5. 

— The scale-height in the solar circle is h z (Ro) — 285l^ 2 
pc. 

— The scale-length of the space density in the plane is 
H = 1.97j£i! kpc. 

— the scale-length of the surface density is ha = 3.3+q'| 
kpc. 

The errors do not include systematic uncertainties, 
which we know to be lower than a 10%. For the flare and 
the warp: 

— The variation of the scale-height due to the flare fol- 

R - R Q 

lows roughly a law h z (R) w h z (R Q )e( 12 - ■» fl i*r c » k P c 
(for R <~ 15 kpc). 

— The warp moves the mid plane of the disc to the height 
z w = 1.2 x IQ- 3 R(kpc) 5 - 25 sin(0 + 5°) pc (for R <~ 13 
kpc), where is the galactocentric angle (</>© = 0). 

The results presented here constitute an important 
step in constraining of the disc structure, especially for 
the stellar populations in outer disc, which were previ- 
ously not very well known. Theoretical implications are 
not discussed here but it is expected that these results 
will constitute a valuable test for many theories about the 
formation and evolution of the galactic discs, as well as 
studies about the formation of warps, flares. 
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Fig. 14. Data and fit of star counts. The dashed line is the best non-warped model {H = 1.91 kpc, h z {R®) = 300 pc, 
hR, flare = 4.60 kpc). The solid line is the best warped model (H — 1.97 kpc, h z (R Q ) — 285 pc, hRji are — 5.00 kpc; 
e w = 5.25, (f> w = -5°, C w = 2.1 x 10" 19 pc). 



